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ABSTRACT 

An analysis of the kinematics of 412 stars at 1-4 kpc from the Galactic mid- 
plane by Moni Bidin et al. (2012b) has claimed to derive a local density of dark 
matter that is an order of magnitude below standard expectations. We show 
that this result is incorrect and that it arises from the assumption that the 
mean azimuthal velocity of the stellar tracers is independent of Galactocentric 
radius at all heights. We substitute the assumption, supported by data, that the 
circular speed is independent of radius in the mid-plane. We demonstrate that the 
assumption of constant mean azimuthal velocity is implausible by showing that it 
requires the circular velocity to drop more steeply than allowed by any plausible 
mass model, with or without dark matter, at large heights above the mid-plane. 
Using the approximation that the circular velocity curve is flat in the mid-plane, 
we find that the data imply a local dark-matter density of 0.008±0.003 M pc~ 3 = 
0.3±0.1 GeV cm -3 , fully consistent with standard estimates of this quantity. This 
is the most robust direct measurement of the local dark-matter density to date. 

Subject headings: Galaxy: disk — Galaxy: fundamental parameters — Galaxy: 
halo — Galaxy: kinematics and dynamics — Galaxy: solar neighborhood - 
Galaxy: structure 



1. Introduction 



The observed flatness of the Milky Way's circular-velocity curve at Galactocentric dis- 
tances larger than 20 kpc (e.g., Xue et al. 2008) shows that the visible Galactic disk is 
embedded in a massive dark halo. The disk is composed of gas and stars (baryons), while 
the dark halo is believed to be composed of non-baryonic matter of unknown nature. Despite 
the dominance of the dark halo in the outer parts of the Milky Way, it remains unclear from 
direct measurements whether there is any need for a substantial amount of dark matter to 
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explain the circular- velocity curve interior to the solar radius, Ro = 8kpc. Quantitatively, 
the fraction of the total radial force at the solar radius that is due to the disk could be 
as high as 90% (Sackett 1997). A promising avenue for constraining the local density of 
dark matter is through a determination of the dependence of the gravitational potential 
on distance above the mid-plane of the disk ("height"), from measuring the kinematics of 
stars (e.g., Kapteyn 1922; Oort 1932; Bahcall 1984). However, a major obstacle is that the 
uncertainty in the amount of baryonic matter in the disk makes it hard to determine the 
relative contributions from dark and baryonic matter to the density near the mid-plane. 
Studies limited to heights < 150 pc are consistent with no dark matter near the Sun, but 
they cannot exclude the amount of dark matter expected for a standard dark-matter halo 
(« 0.01 M pc" 3 = 0.38 GeV cm" 3 ; Creze et al. 1998; Holmberg & Flynn 2000). 

The contributions from baryonic and dark matter can be disentangled by measuring the 
gravitational potential out to larger heights. At heights of several times the disk thickness, 
the dark halo and the baryonic disk contributions to the potential have a different vertical 
dependence (e.g., Kuijken & Gilmore 1989; Garbari et al. 2011). In particular, most of the 
disk mass lies below ~ 1.5 kpc, so above this height the disk contribution to the integrated 
surface density = J_ z dz p(z) is roughly constant with height and the disk potential 

varies as $ oc \Z\, while the dark halo contributes oc \Z\ and $(Z) oc Z 2 . Thus, any 

measured increase in the surface density at \Z\ > 2 kpc must be due to the dark halo, and 
the expected increase is ~ 20 M & pc -2 for each kpc (for the standard value of the local dark 
matter density of 0.01 M pc -3 ; see above). Determinations of the surface density at ~ 1 kpc 
from the plane typically find values of E(lkpc) = 70 to 8OM pc -2 (Kuijken & Gilmore 
1991; Siebert et al. 2003; Holmberg & Flynn 2004) while the baryonic contribution is esti- 
mated to be around 50 to 6OM pc -2 (Holmberg & Flynn 2004; Bovy et al. 2012a). Thus, 
these studies are consistent with the expected dark-matter density, although they do not go 
to large enough heights to detect the dark matter unambiguously. 

A recent study of the kinematics of stars with heights 1.5 < \Z\ < 4 kpc by 
Moni Bidin et al. (2012b, MB12 hereafter) claims to have measured the surface den- 
sity at these heights in a model-independent manner and found it to be constant, such that 
the local density of dark matter must be smaller than lO~ 3 M pc~ 3 (< 0.04 GeV cm -3 ). If 
true, this would imply that dark matter is less abundant in the solar neighborhood than 
expected, by at least an order of magnitude. This would also shift experimental limits on the 
cross-section for elastic scattering between dark-matter particles and baryons by an order 
of magnitude, although this is far less than the uncertainty in the predicted cross-section. 

In this paper, we show that the analysis used by MB 12 is flawed. The main error is that 
they assume that the mean azimuthal (or rotational) velocity V of their tracer population 
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is independent of Galactocentric cylindrical radius R at all heights (i.e., V(R,Z) = V(Z)). 
This assumption is not supported by the data, which instead imply only that the circular 
speed V c is independent of radius in the mid-plane (e.g., Gunn et al. 1979; Feast & Whitelock 
1997). In the solar neighborhood, the circular speed is larger by > 35 km s _1 than the mean 
azimuthal velocity for the warm tracer population used by MB12, a phenomenon known as 
asymmetric drift. The asymmetric drift is expected to vary with R — although this variation 
cannot be measured for the sample of MB12 as the data do not span a large enough range 
in R — so the assumptions that V and V c are independent of radius are not compatible. 
In the absence of a measurement of the (R, Z) dependence of V for the tracers, we show 
that the assumption of an i?-independent V at all heights Z is highly implausible. By using 
instead the data-driven assumption that the circular speed is independent of radius at Z — 0, 
we demonstrate that the measurements and analysis of MB 12 are fully consistent with the 
standard estimate of the local dark matter density, approximately 0.01 M pc -3 , and indeed 
provide the best available direct measurement of this quantity. 

The outline of this paper is as follows. In § 2 we discuss the assumption made by MB 12 
that the mean azimuthal velocity of the tracer population is independent of R. In § 3, we 
derive the surface density as a function of Z using the assumption that the circular-velocity 
curve is flat in the mid-plane and we discuss the effect of relaxing this approximation. In 
§ 4, we calculate the surface density using the data of MB12 but using the approximation 
of constant circular speed and show that this leads to a surface density at 1.5 < \Z\ < 4kpc 
that increases with height at the rate expected for a dark halo. Our conclusions are in § 5. 

We follow MB12 in defining U, V, and W as the radial, azimuthal, and vertical velocities 
in cylindrical coordinates and the inertial Galactocentric reference frame. Note that this 
definition is non-standard, since (i) U, V and W are normally defined with respect to the 
Local Standard of Rest; (ii) U is typically the velocity toward the Galactic center; (iii) our 
conventions require either that the coordinate system is left-handed or that the positive 
z-axis points to the South Galactic Pole. However, for the purpose of this paper these 
distinctions do not matter. We use Ro = 8kpc, as do MB12, and local circular speed 
V c = 220 km s" 1 throughout this paper, assumptions which are consistent with the latest 
measurements (Bovy et al. 2009). We abbreviate the surface density at the solar radius as 
Z(Z)=Z(Ro,Z) = f* z dzp(Ro,z). 

2. Mean azimuthal velocity versus circular velocity 

The analysis of MB12 uses the Poisson equation to calculate the surface density as a 
function of height Z at the solar radius Rq using the radial force Fr and the vertical force 
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Fz (eq. 18). To estimate Fr, MB12 make use of the radial Jeans equation 

d*(R,Z) 1 d(ua 2 ) 1 d{va 2 w ) a* - a\ - V 2 
F * R > Z > = dR— = u^R- + u^Z— + R ' (1) 

where <3> is the gravitational potential, v is the tracer-density profile, afj and a v are the radial 
and azimuthal velocity dispersions squared, <Ju W is the off-diagonal radial-vertical entry of 
the dispersion-squared matrix, and V is the mean azimuthal velocity; all of these quantities 
are functions of R and Z. We have assumed a steady state such that time derivatives 
vanish and the mean radial motion U is zero. In simplifying equation (1), MB12 use their 
assumption VIII: "The rotation curve is locally flat in the volume under study", which MB12 
express as 

Moni Bidin et al. (2010) made the same assumption, although this was not explicitly stated. 
We show in this section that this is an unreasonable assumption. 

The mean azimuthal velocity of a population of stars differs from the circular velocity 
due to the asymmetric drift. This offset arises because both the density of stars and the 
velocity dispersion typically decline with radius. This means that more stars with guiding 
centers at R < Rq are passing through the solar neighborhood than stars with guiding centers 
R > Rq] the former are on the outer parts of their orbits, where their azimuthal velocity 
is less than the circular velocity. Equation (1) is typically used to estimate the asymmetric 
drift (Binney & Tremaine 2008): since V 2 = Rd^/dR we find 



V 2 -V 2 = a 2 v -a 2 u 



i din (yap 1 R d [ya 2 uw ) 



dlnR v al dZ 



(3) 



Now assume that (i) the dispersions-squared (cry, cr v , cr1j W ) decline exponentially with 
radius with scale length h a ; (ii) the tracer density is an exponential function of radius and 
height with scale lengths h R and h z , that is, v(R, Z) oc exp(— R/Hr — \Z\jhz). The second 
assumption is only accurate at heights above a few hundred pc; closer to the mid-plane the 
exponential form is not accurate but this is not a concern since we are interested in the region 
\Z\ > 1 kpc. These assumptions were also made by MB12 and in addition they assumed that 
h a = hR. Equation (3) then becomes 



V c 2 -V 2 = a 2 v + a 2 u 



^ ' h R + h t ' 



+ Y z auw - R ^z~- (4) 



To evaluate this, we use the expressions for the moments of the velocity distribution from 
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Moni Bidin et al. (2012a) 



au(Ro, Z) = (82.9 ± 3.2) + (6.3 ± 1.1) ■ (|Z|/kpc - 2.5) km s" 1 (5) 
a v (Ro, Z) = (62.2 ± 3.1) + (4.1 ± 1.0) ■ (\Z\/kpc - 2.5) km s" 1 (6) 
a w (R , Z) = (40.6 ± 0.8) + (2.7 ± 0.3) • (\Z\/kpc - 2.5) km s" 1 (7) 

and that for (rfj W from MB 12 

a 2 uw (R , Z) = (1522 ± 100) + (366 ± 30) ■ (\Z\/kpc - 2.5) km 2 s" 2 . (8) 

We also take h z = 0.9 kpc and = h a = 3.8 kpc as in MB12, although our own work 
suggests somewhat smaller values Jir = 2 kpc and h a = 3.5 kpc (Bovy et al. 2012c). We find 
that the solution to equation (4) can be fit by the formula 

V 2 - V 2 = (191 km s" 1 ) 2 [1 + 0.19 (\Z\/kpc - 2.5)] , (9) 

with an rms error of less than 2% for \Z\ = 1-4 kpc. Equation (9) implies V c — V(1.5 kpc) = 
80 km s _1 and V c — ^(3.5 kpc) = 150 km s _1 , in good agreement with the corresponding 
measured values in Moni Bidin et al. (2012a), 70 ± 13 km s -1 and 130 ± 16 km s^ 1 . 

We can now ask what equation (4) predicts for the radial variation of V. Taking the 
radial derivative of equation (4), we find that 

r SVc 9V dV (V c 2 - V 2 ) 2 / 1 1\ 1 2 da 2 uw 
"dR- 2V dR= — + au \h- R + K) + ~ ^Z~ 



(K 2 -^)(^-T^)+^-^- (10) 

We can then estimate dV /dR for a flat circular-speed curve, by setting dV c /dR = 0: 
8V 

V— = 110 km s" 1 x 21km s" 1 kpc" 1 [1 + 0.2 (|Z|/kpc - 2.5)] , (11) 
oR 

with an rms error of less than 1% for \Z\ = 1-4 kpc. In this equation, 110 km s _1 and 
21 km s" 1 kpc" 1 are V and dV/dR at \Z\ = 2.5 kpc. We find that dV/dR is 7km s" 1 kpc" 1 
at Z = 0, growing to 11km s _1 kpc -1 at \Z\ = lkpc and 40 km s _1 kpc" 1 at \Z\ = 3.5 kpc. 
If we use hfj = 3.5 kpc (Bovy et al. 2012c), the gradients are larger by about 20%. 

For comparison, MB12 estimate that a gradient dV /dR = 17 km s" 1 kpc" 1 is needed 
to make their analysis consistent with the expected amount of dark matter, which is close 
to our estimate of 21km s~ 1 kpc~ 1 . MB12 dismiss this possibility, apparently because they 
confuse constraints on dV c /dR with constraints on dV/dR, as there are, in fact, no direct 
observational constraints on dV /dR. 
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While dV/dR is not measured for the MB12 sample, constraints on this gradient do 
exist for a slightly more metal-rich sample of stars at moderate distances from the mid- 
plane, which we can use to show that our estimate of dV/dR using equation (10) agrees 
with observations. Casetti-Dinescu et al. (2011) study a sample of metal-rich red clump 
stars at 1.0 kpc < \Z\ < 2.5 kpc, selected using 0.6 < J — K s < 0.7, while the MB12 
data were selected using 0.7 < J — K s < 1.1 to obtain a lower- metallicity, larger- distance 
sample. Casetti-Dinescu et al. (2011) report that (au, &v, V) = (60, 42, 180) km s -1 at \Z\ « 
lkpc, such that equation (10) assuming a flat circular-speed curve predicts that dV/dR = 
6.5 km s" 1 kpc" 1 ; similarly, they find (au, ay, V) = (80, 60, 150) km s" 1 at \Z\ m 1.5 kpc 
and (au,(Ty,V) = (90, 65, 140) km s" 1 at \Z\ « 2kpc, such that we predict that dV/dR = 
13 km s _1 kpc _1 and 15 km s _1 kpc, respectively. These predictions are in good agreement 
with the measurements of Casetti-Dinescu et al. (2011), who find 6.0 ± 1.5km s -1 kpc -1 , 
13.0 ±4.5 km s _1 kpc~\ and 12.0 ±8.1 km s" 1 kpc" 1 at these \Z\. 

The discussion so far in this section has assumed that V c is constant with Z and R, 
whereas observations only show that V c is independent of R at Z = 0. The circular velocity 
declines with \Z\ for any reasonable local mass distribution, as the following argument shows. 
As V^? = —RF R , we can write 

dV 2 d dF 7 

where we have used the fact that Fr and Fz are both derivatives of the potential $, such 
that 

dF R d_(d±\ d (d®\ = 8F Z 

dZ dZ\dRJ dR\dZJ OR' 1 j 

As shown below in § 3, to a good approximation Fz = — 2ttGT,(R, Z) (for Z > 0), where 
E(i?, Z) is the surface density. Therefore 

2 V c § . 2nGR ^§^1 = -2nG^- E (j R, Z) , (14) 
oZ oR /is 

assuming that H(R,Z) declines exponentially with radius with scale length h^, which we 
assume is equal to 3.5 kpc (Bovy et al. 2012c). This means that near the mid-plane, where 
the density is approximately 0.1M Q pc~ 3 , 

dV r . 



-2.8 km s" 1 kpc" 1 ( ) \Z\ < h z 



oz 

while for \Z\ > 1.5 kpc, where E(R , Z) ~ 50 M pc~ 2 



(15) 



8Vr 



dZ 



c 



-7km s" 1 kpc" 1 \Z\ > h z . (16) 
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Taken together, these results imply that V c likely does not decrease by more than 30 km s 
out to \Z\ = 4kpc. Such a decrease does not change the conclusions of this section. 

A different way of illustrating the inconsistency of MB12's assumptions is to ask what 
the radial behavior of V c has to be such that dV/dR = 0. The necessary dV c /dR decreases 
from approximately —5 km s _1 kpc -1 at Z — 0, which is consistent with local measurements 
of the slope of the circular-speed curve, to —14 km s _1 kpc -1 at \Z\ = 4 kpc; for h a = 3.5 kpc, 
the gradients are larger in absolute value by 20%. Such a steep drop of the circular velocity 
with R is about the value of a Keplerian drop-off, —^V c /R = —14 km s" 1 kpc^ 1 . Such steep 
gradients are inconsistent with observational evidence: (i) in the limiting case where there 
is no dark-matter halo and all the mass is in an exponential disk (with parameters given in 
Figure 1), \dV c /dR\ is only 3km s~ 1 kpc~ 1 at \Z\ = 4 kpc, and if a dark halo is present the 
value is even smaller; (ii) we show below in § 3 that dV c /dR increases as \Z\ grows, whereas 
maintaining constant V requires that dV c /dR decreases with \Z\. 



3. Poisson equation at large heights 

Having shown that the assumption of a radially constant V is suspect, we now ask what 
the MB12 data do have to say about the surface density at large height \Z\. Following MB12, 
we start from the Poisson equation in cylindrical coordinates, 



£(#, Z) 



2uG 



1 d(RF R ) 
dZ R^R- + Fz{R > Z) 



(17) 



where we have assumed symmetry around the Z = plane. Fz(R, Z) is obtained from the 
steady-state vertical Jeans equation 

d$(R,Z) 1 d{Rval w ) 

Fz(i? ' Z) = dZ— = u^Z- + R^ dR ■ (18) 

One can then proceed (e.g., Kuijken & Gilmore 1989) by approximating the integrand in 
the first term in square brackets in equation (17) by its value in the plane. This is zero 
for a flat circular-speed curve, as RF R = —V^. Kuijken & Gilmore (1989) show that this 
approximation is good to a few percent at \Z\ < 1.5 kpc for any plausible mass distribution. 
We now revisit this question at larger heights. 

Figure 1 shows the error introduced by approximating d (RFr) / dR by its value at 
Z — 0, for three different mass distributions: an exponential disk with parameters that 
are representative of the Milky Way's disk (Bovy et al. 2012a,b), a Navarro- Frenk- White 
(NFW; Navarro et al. 1997) halo, and a combination of the two in which the disk provides 
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85 % of the radial force at Rq (this yields a flat circular-speed curve at Rq). This figure shows 
the fractional difference between the true surface density and that calculated from equation 
(17) by approximating d (RF R ) /OR by its value in the plane. The errors introduced by 
this approximation can be as large as 15 % at 4 kpc above the plane. Thus a systematic 
uncertainty of 10 to 20% in E(4 kpc) is unavoidable without fully modeling the potential. 

It is, however, possible to show that neglecting the Z-dependence of dV c /dR almost 
always leads to an underestimate of the surface density, so this approximation gives a robust 
lower limit on the surface density. We start by writing the Z derivative of the integrand in 
equation (17) as 

d (l d{RF R )\ 8 (F R , dF R \ 



dZ \R OR J dZ \R OR 
If we again use the fact that dF R /dZ = 8F z /dR, we find that 



8 ( 1 d(RF R )\ 1 dF z d 2 F z 



dZ \R dR J R OR dR 2 ' 

If the first term in square brackets in equation (17) is neglected, we have F z = — 27rGS(i?, Z) 
(at Z > 0). This approximation is quite accurate if there is substantial mass in a thin disk 
and \Z\ is small, or if the circular-speed curve is nearly flat, since then RF R is almost 
independent of radius (see Figure 1). If we then drop this term and approximate the radial 
derivatives then as those of the exponential disk, equation (20) becomes 



d (\ d(RF R ) \ _ 2txG f R 
dZ\R dR J ~ Rh s V ~ 



E(R,Z). (21) 



Since Ro/h^ is 2 or more in any reasonable model for the Milky Way, the right side of this 
equation is negative, which in turn implies that the integrand in equation (17) decreases with 
increasing \Z\. Thus replacing F R by its value in the mid-plane leads to an underestimate 
of the surface density, and dropping the integral in equation (17) when the circular-speed 
curve is locally flat also leads to an underestimate of the surface density. 

We can further work out equation (21) to find the change in the radial slope of the 
circular-speed curve with height 

8 ravA nG fRA^^ (22) 



dZ \dRJ hxV c \h s 

a result that can also be derived by differentiating equation (14). Using the same reasoning 
that led to equations (15-16), we find 



° 0.45 km s" 1 kpc" 2 ( 

dZ\dRJ 1 ViOOpc 



\Z\<h z , (23) 
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and 

— (-^j ^l.lkms^kpc" 2 , \Z\~>h z . (24) 

Thus, the radial gradient of the circular speed remains close to its value in the plane through- 
out the region \Z\ < 4kpc that we are investigating. 



4. Moni Bidin et al. (2012) revisited 

We have shown that approximating the integrand in equation (17) as constant with Z 
leads to a robust lower limit on the surface density as a function of Z. Armed with this result, 
we can revisit the MB 12 analysis. Since the circular-speed curve is flat in the mid-plane, 
the integrand is zero in the mid-plane, so the approximation that the integrand is constant 
implies that the integral can be neglected. Substituting equation (18) into equation (17) we 
find 



E(Z) 



1 



2ttG 



daL . / 1 1 1 



2 , ™W , 2 



(25) 



h z w dZ uw \R h R h 
where we distinguish between the radial scale length Hr of the tracer population and the 
radial scale length h a of the dispersion-squared. MB12 assumed h R = h„ = 3.8 kpc, but 
recent direct measurements of the thicker populations of stars in the Milky Way have shown 
that their radial scale length is much shorter than that of the dominant, thinner components, 
such that Iir = 2 kpc is more accurate; we also prefer a shorter scale length for the dispersions, 
h a = 3.5 kpc (Bovy et al. 2012b, c). We will nevertheless show results mostly for the values 
h.R = h a = 3.8 kpc assumed by MB12. We also follow MB12 in assuming hz = 900 pc (see 
below). Note that the surface density as calculated using equation (25) does not depend on 
the value of the local circular speed V c , and our estimate of the local dark matter density is 
therefore not affected by the uncertainty in V c . 

Using the kinematics from Moni Bidin et al. (2012a) and MB12 in equation (25), we 
get the curve labeled as "correct dV c /dR = approximation" in Figure 2. As we have 
discussed, this gives a lower limit to the surface density. The gray shaded region in the 
Figure shows the range of surface densities obtained after including the radial integral in 
equation (17), using the bottom and top curves in Figure 1. These can be compared to the 
curve labeled "incorrect dV/dR = approximation", which is the curve presented in MB12 
as their primary result (their Figure 1). 

We can also compare these curves with Galactic mass models in the literature, showing 
the same mass models as MB12. The estimated contribution to from baryonic matter 
(VIS for "visual" in Figure 2) is composed of a thin 13 M Q pc~ 2 layer of interstellar medium, 
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plus a stellar halo and "thick" and "thin" disk components with parameters taken from 
Juric et al. (2008); in this model the normalization is chosen so that Sdisk(l-1 kpc) = 4OM 
pc" 2 (Holmberg & Flynn 2004). Figure 2 also shows the expected surface density when 
various dark-matter halo models are added to this baryonic mass. The OM halo model has 
a profile 



with R c = 8.01 kpc and p c = 0.0103 M pc" 3 (Oiling & Merrifield 2001), such that p 0jDM = 
0.0084 M & pc" 3 . The other halo models were all taken from Weber & de Boer (2010) and 
have densities following 



where r is the Galactocentric spherical radius. The Standard Halo Model (SHM) has an 
NFW profile (a = P = 1, 7 = 2) with R c = 10.8 kpc and p 0>DM = 0.0084 M pc" 3 . The 
N97 model also has an NFW profile with R c = 20 kpc and p 0>D M = 0.0061 M pc" 3 . The 
MIN model has a profile as in equation (26) (i.e., a — 0, /? = 2, 7 = 1) with R c = 5 kpc and 
p c = 0.019 M Q pc" 3 (Weber & de Boer 2010), such that p 0jDM = 0.0053 M pc" 3 . All of 
these models were constrained by assuming that V C (R ) = 244 ± 10 km s~\ except for the 
OM model which has V c (Ro) = 220 km s" 1 . The difference between the predictions of these 
models for £(Z) in Figure 2 is mainly due to their different value for p ,dm- 

It is clear from Figure 2 that the MB12 data are fully consistent with the predictions from 
several standard dark matter models when using the correct assumption dV c /dR = 0. In 
particular, the values of the surface density throughout the range \Z\ — 1-4 kpc are consistent 
with the standard halo model (SHM) which has pdm — 0.01 M pc" 3 = 0.38 GeV cm" 3 . 
Figure 2 also shows the effect of using the more appropriate radial scale length Hr \ — 2 kpc, 
which increases the surface density even further. The slope of the measured Ti(Z) implies a 
minimum dark-matter density of p DM = 0.007±0.001 M pc" 3 = 0.27±0.04 GeV cm" 3 for 
h R = 3.8 kpc, and p DM = 0.0085±0.0015 M pc" 3 = 0.32±0.06 GeV cm" 3 for h R = 2.0 kpc. 
These uncertainties are statistical and do not include the systematic uncertainty associated 
with the approximation that the circular-velocity curve is flat at all heights, which adds 
at most 0.1 GeV cm" 3 (see Figure 1), or that associated with the measurement of the 
dispersions (see, e.g., the erratic behavior of afj W in figure 8 of Moni Bidin et al. 2012a), 
or the assumption that the tracers follow an exponential distribution in radius and height 
with the assumed scale lengths. These systematic uncertainties are at least of the same 
magnitude as the statistical uncertainties. For example, if instead of hz = 900 pc we use 
hz = 700 pc, which is the mass-weighted mean scale height in the range 2 to 4 kpc using 
the measurements of Bovy et al. (2012a,b), pdm increases by 0.001 M pc" 3 ; using a varying 




(26) 




(27) 
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mass- weighted hz as a function of \Z\ as given by those same measurements changes pdm 
by 0.002 M Q pc -3 . Taking these systematics into account, our best estimate for the local 
dark-halo density is p DM = 0.008 ± 0.003 M & pc" 3 = 0.3 ± 0.1 GeV cm" 3 . 

5. Conclusions 

In this paper, we have shown that the assumption of a radially constant mean azimuthal 
velocity for the stellar tracers used by MB12 is physically implausible. This assumption is 
the reason why the MB12 analysis finds a constant at 2 < \Z\ < 4 kpc, in seeming 

contradiction with the standard expectation of a dark-matter density pdm ~ 0.01 M pc~ 3 = 
0.38 GeV cm -3 in this region. The mean azimuthal velocity V of the warm tracer population 
of stars used by MB12 — stars that are part of the thicker components of the Milky Way disk 
that reach large heights above and below the plane — is significantly different from the circular 
velocity V c , with V < 0.5^ at \Z\ = 3 kpc. We have shown in § 2 that for a circular-speed 
curve that is close to flat, the expected dV/dR reaches tens of km s _1 kpc _1 at a few kpc 
from the Galactic mid-plane. Indeed, the assumption that dV/dR = requires that the 
circular-speed curve falls off in a Keplerian manner at a few kpc above the plane, in clear 
contradiction with observations. 

We derived an alternative formula for that assumes that the circular-velocity 

curve is flat at all heights above the plane and we showed that this approximation leads to 
a lower limit for all plausible mass distributions. This approximation sidesteps the issue of 
the unknown radial trend of V and as such makes no assumptions about it. Applying this 
alternative formula, we find that the MB 12 data give a lower limit that is fully consistent with 
the standard local density of dark matter of p^M ~ 0.01 M pc~ 3 , and that they imply a local 
dark-matter density of 0.008 ± 0.003 M Q pc -3 = 0.3 ± 0.1 GeV cm -3 , where the error bars 
include both statistical and less well-known systematic errors. Therefore, our analysis shows 
that the locally measured density of dark matter is consistent with that extrapolated from 
halo models constrained at Galactocentric distances r > 20 kpc: for example, extrapolating 
the best-fit halo profiles of Xue et al. (2008) and Deason et al. (2012), obtained from fitting 
halo stars with 20 < r < 60 kpc, gives pdm ~ 0.006 to 0.011 M Q pc~ 3 . 

The breakdown of the assumptions made in this simple, "model-independent" Jeans 
analysis are such that the measurement has a systematic uncertainty reaching 10 to 20% 
at \Z\ = 4 kpc. Therefore, a precise determination of the local dark matter density from 
observations at large Z using the Jeans analysis of MB12 requires data that span a wide range 
in R such that the radial gradient of the velocity moments, foremost V, can be determined. 
The Gaia mission (Perryman et al. 2001) will provide such measurements in the near future. 
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Currently, data that span multiple kpc in R and Z are available from the SDSS/SEGUE 
survey (Yanny et al. 2009). In contrast to the data from Moni Bidin et al. (2012a), the 
SEGUE volume selection is known in detail (e.g., Bovy et al. 2012b) so both the spatial 
structure and kinematics can be obtained for the same set of tracer stars. The SEGUE data 
also have the advantage that they contain information on the elemental abundances of each 
star, allowing the tracers to be divided into sets of stars with simple distribution functions 
(e.g., Bovy et al. 2012b,c). 

Finally, we note that our estimate is of the mean halo density between 1 and 4 kpc 
above the Galactic mid-plane. The halo density at the Sun, which is the relevant quantity 
for direct dark-matter detection experiments, is likely to be larger because of two effects. 
The density in the mid-plane for a spherical NFW halo with a scale radius of 22.25 kpc 
(Xue et al. 2008) is 7% larger than at \Z\ = 2.5 kpc. Besides this purely geometric effect, the 
gravitational influence of the disk further increases the mid-plane dark-matter density. An 
isothermal halo with isotropic velocity dispersion a has a density p oc f dv e~ E l a2 = e - */"" 2 , 
so we expect that 

Pdm(\Z\) 27rG£ disk |Z| _ _ Q 2Q E disk / 130 km s -1 \ 2 \Z\ 

p DM (0) ' a 2 ' 5OM pc- 2 \ a ) 2.5 kpc' [ ' 

These two effects imply that the dark-matter density in the mid-plane is enhanced over the 
value derived in this paper by about 30%. This agrees with the estimated enhancement in 
an iV-body simulation by Garbari et al. (2011). 

It is a pleasure to thank the anonymous referee, Dana Casetti-Dinescu, Christian Moni 
Bidin, Jerry Ostriker, Justin Read, Hans- Walter Rix, Martin Smith, and Nadia Zakamska 
for helpful comments. Support for Program number HST-HF-51285.01-A was provided by 
NASA through a Hubble Fellowship grant from the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, Incorporated, 
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Fig. 1. — Fractional difference between the true surface density S true (i? , Z) and that 
obtained by approximating the integrand in equation (17) by its value in the plane, 
E approx (i?o, Z). Shown are an exponential disk, a spherical NFW halo, and a combination of 
the two that has a circular-speed curve that is flat near Rq = 8 kpc. 
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Fig. 2. — The surface density as a function of height using the invalid assumption that 
dV/dR = (lower black curve) and the more realistic assumption that dV c /dR = in the 
mid-plane (upper black curve). The latter assumption is shown in § 3 to give a robust lower 
limit to the surface density. The dashed curve shows the effect of reducing the radial scale 
length of the tracer from MB12's value Hr = 3.8 kpc to the more likely value of 2 kpc. Also 
shown as the gray band is the range of surface densities that results from applying the lower 
and upper curves in Figure 1 to correct the approximation that V c is independent of height; 
a similar gray band would apply to the dashed curve. 68 % uncertainty intervals on the 
observed surface density are shown at a few representative points. The curves representing 
estimates of the visible matter ('VIS') and the predictions of various dark-matter halo models 
('OM','SHM','N97', and 'MIN'), defined in § 4, are the same as in Figure 1 of MB12. 



